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Abstract 



If a heavy Dirac monopole exists, the light to light scattering below the 
monopole production threshold is enhanced due to the strong coupling of 
monopoles to photons. This effect could be observable in the collision of 
virtual photons at proton colliders. At the Tevatron it will be seen as pair 
production of photons with energies 200-400 GeV and roughly compensated 
transverse momenta 100-400 GeV/c. This effect could be seen at monopole 
masses about 1-2.5 TeV at the upgraded Tevatron and 7.4-19 TeV at LHC 
depending on monopole spin. 
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I. INTRODUCTION 



The magnetic charge (monopole) was introduced into particle theory by Dirac P| (see 
also [0) to restore symmetry among electricity and magnetism. The idea of a monopole 
is very attractive in order to explain the mysterious quantization of the electric charge of 
particles. Therefore, the search for a monopole is essential despite the fact that there is no 
place for it in the standard description of our world. 

The Dirac-Schwinger monopoles, discussed here, are point-like particles. They differ 
strongly from non-local monopoles in the context of gauge theories, first considered by 
Polyakov and 't Hooft || . Below we assume that a few monopoles exist in the Universe and 
they were not yet observed due to their high mass. 

Our basic idea is simple: The existence of monopoles provides for a 77 elastic scattering 
at large angles, that is sufficiently strong below the monopole production threshold. This 
effect is observable at colliders with energies smaller than monopole masses. Such a method 
for discovering monopoles has been suggested in Refs. A similar idea for the process 

e + e" — ► Z — > 37 has been considered recently and tested at LEP ||. Our paper deals 
with detailed calculations of this effect at hadron colliders. Throughout the paper we denote 
the monopole mass by M and its spin by Jm (we assume a definite spin of the monopole), 
lo is the characteristic photon energy (typically - the 77 c.m.s. energy). 

A theory with two point-like charges, electric and magnetic, cannot be standard QEDQ. 
According to Refs. |I]fJ] (see also the electromagnetic field in such a theory is described 
by a vector potential having the Dirac string or some of its surrogate. To have unambiguous 
results, the elementary electric and magnetic charges e and g ought to be quantified so that| 

2™ . , 

9 = , n = ±l, ±2,... 1 

e 

with a = e 2 /(47r) = 1/137 =>- a g = g 2 /(A7c) = n 2 /(4a). Unfortunately, the explicit form of 
such a theory is still not yet known. 

We are interested in the energy region below the monopole production threshold. Here 
the interactions of photons via virtual monopoles are considered as main effect. The effective 
expansion parameter is of the order of 

9eS = r— = < 1- (2) 

VA^M 2^M 1 ; 



In this region general considerations like gauge invariance, threshold behavior, etc. together 
with a perturbative approach allow us to believe that a perturbation theory analogous to 
standard QED can be applied (certainly, in lowest nontrivial order only). In that case 
QED-like calculations should be valid, both in tree approximation and at one loop level. 
For one loop we can assume a Wick rotation into Euclidean region. As it is well known, the 



1 According to ||, this theory even could violate U(l) local gauge invariance of standard QED. 
We are thankful to Dr. He for clarification of this point. 

2 In the Schwinger theory [0] the quantity n should be even, n = ±2, ±4, . . . 
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integration over the loop momentum p is convergent due to gauge invariance. Therefore, the 
integration region is limited by virtualities p 2 ~ u 2 , where the effective expansion parameter 
is small, and QED results are valid. 

The estimate of the effective coupling constant is supported by estimating qualitatively 
the cross section of 77 — > 77 scattering via a monopole loop (Fig. [I]) in the 77 c.m.s. 
(without fixing the form of the interaction). Due to gauge invariance, each photon leg of the 
matrix element yields a factor u. On dimensional grounds this factor has to be divided by the 
monopole mass M. Additionally, the magnetic charge g has to be associated to each vertex. 
Therefore, the amplitude M. oc (goo/M) 4 and the cross section isQ a oc {1/uj 2 ) (gu/M) 8 . 

The corresponding effective Lagrangian (of Heisenberg-Euler type) is 



cff 



+ 



36(4tt) 2 M 4 



(3 + +(3- 



P+-P- 



+ ... (3) 



with the electromagnetic field strength tensor F^ v and F^ v = e^ va ^ F a ^j2. An additional 
gauge fixing term has to be introduced to invert the photon propagator. The constants 
At depend on the monopole spin (a numerical coefficient is introduced to simplify the final 
expressions). Based on arguments mentioned earlier (related to Wick rotation) we use here 
the coefficients f3± obtained in QED. Their values have been found for different values of 
the monopole spin J M in (J M = 1/2), (TJ (J M = 0), [l|Jl| (J M = 1). We use their 
combination P = j3\ + 2f3 2 _. 

Let us denote the standard Mandelstam variables for the 77 scattering by s, t, u with 
s +t + u = (sis the effective mass of produced 77 system squared and t varies within the 
interval ( 0, — s/2)). Furthermore, 6 is the scattering angle in c.m.s. of the photons with 
< 9 < 7r/2. In our case s = 4tu 2 , i = -2u 2 (l - cosQ). 

The cross section for 77 scattering via a monopole loop at s M 2 is given by the 
expressions 



a 




R 



11 



28P 



405vr \2 v / aM / 



(4) 



da 



5(3 



cos 



56 



-a tot dcosfl 



5 fs 2 + i 2 + u 2 \ 2 



dt 



a 



(5) 



f 0.085, J M = 

P = (3 2 + + 2f3 2 _ = \ 1.39, J M = 1/2 . (6) 
{ 159, J M = 1 

Note that the result strongly depends on the monopole spin Jm- 



3 If some monopole - antimonopole bound state with mass <C 2M exists, it is wide enough (since 
g 3> 1), and it gives additional more strong light to light scattering. We don't speculate about this 
opportunity in detail. 
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These expressions are correct up to contributions of order 0(g 2 s ). We expect that with 
growing g eS the increase of the cross section as a function of energy stops and the production 
of a larger (even) number of photons becomes essential. The effective parameter here is 
expected to be less than g 2 s . For example, in QED with spin 1/2 fermions, the ratio of 
coefficients of the 67 to 47 operators in the effective Lagrangian for 77 interactions is about 
0.3g 2 s BUI. Besides, one can expect that below the threshold (a; < M) the expressions (|J) 
correctly describe the sum of cross sections of processes 77 — * 27, 47, 67, . . . with increasing 
multiple photon production at higher g e g. 



II. 77 PRODUCTION VIA A MONOPOLE LOOP AT HADRON COLLIDERS 

The cross section for the production of two photons via a virtual heavy monopole loop in 
pp or pp collisions is a convolution of the 77 cross section (|^) with the photon fluxes arising 
from the colliding protons 



a 



pp—*'yyX 



- n 1 {u) ll Q l )n 2 {uj2,Q 2 ) 2 2 rfff 77 ^ 77 . (7) 

7T/ J UJ\ LU 2 Qi Q 2 



Since in the dominant integration region Q 2 <C s, the transverse motion of initial photons 
can be neglected with good accuracy 

The relevant scale for the virtuality dependence of the 77 — > 77 subprocess cross section 
is given by the unique inner parameter of monopole loop - the monopole mass. Therefore, 
this dependence appears in the form of the quantity Q 2 /M 2 . Since Q 2 <C s <C M 2 , it is 
safely neglected below. 

Taking these approximations into account the cross section (0) can be written in factor- 
ized form (valid for Q 2 <C s). The photon flux densities rij depend on their "own" photon 
variables only 



rii =ni{yi= Q 2 ) , s = 4uiU 2 (8) 



and we have for the cross section 

'a\ 2 r dyi , ,dy 2 I( N _,_ _ / <\ x 2 
^) J yi 1 V2 



a 



pp—fyyX 



I ^m^m^y^n - i-) 2 RE'N 2 {E) , (9) 
J y\ y 2 Vvr/ 

N(E) = J y 2 f{y)dy , f(y) = J n(y, Q 2 )^ . (10) 

The photon flux density arising from one proton is a sum over elastic and inelastic 
contributions: 

n(y, Q 2 ) = (D cl + An) + j(C cl + C- m ) . (11) 

In the individual contributions the quantity Q 2 is limited kinematically from below (m p is 
the proton mass): 

Q 2 > QLn = {M 2 X -m 2 p ) T ^ + m 2 p ^— . (12) 
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Mx is the effective mass of the system produced in the virtual 7*p collision, Mx = m p in 
the elastic case. 

The elastic contribution is written via standard proton form factors Ge and Gm 
C c i(Q 2 ) = G 2 M (Q 2 ) , 

n(n 2, n 4m 2 p G 2 E (Q 2 )+Q 2 G 2 M (Q 2 ) ( Q 2 min \ 

The integral over Q 2 in ([H]) rapidly converges at the upper boundary due to the form factors. 
The integral is saturated at Q 2 values given by the form factor scale. This scale is much 
lower than the other parameters of the problem. Therefore, the upper integration limit can 
be extended to infinity and the elastic contribution to N(E) becomes energy independent. 

The inelastic contribution is written via the proton structure functions F\ , F 2 . C- in and 
-Dj n are integrals over the squared effective mass of the produced system M x . 

We basically start from Eq. (D.4) and Table 8 of Ref. |14j 

C in (Q 2 ) = -^JdM 2 x F 1 (M 2 x ,Q 2 ), 

D hl (Q 2 ) = ^Zfi- I dM 2 x xF 2 (M 2 x ,Q 2 ) ( 1 - . (14) 



Q J -A J -V / ^ Q2 

In order to use the standard representations of the structure functions we change the integra- 
tion variable from M\ to the Bjorken variable x using the relation M\ = m 2 + Q 2 (l — x)/x. 
Next, from inequality dl2|) the lower limit in x is found 

^min = Z 2 2 ;^i2 ' (1^) 

1 - y m p /Q 



yielding the useful relation 



( l -yM 1 -%r) = y(-^-± 



The upper x value is reached at the minimal value for the quantity L = M x — m 2 ~ 2m p m 7r . 
With these transformations we present C m (Q 2 ) and Di n (Q 2 ) in the form 

^max i 2? max , 

C in (Q 2 ) = 2 / -jFi(x, Q 2 ) , D in (Q 2 ) = y I ^(—-i)f 2 (x,Q 2 ). (17) 

^min ^min 

Concerning the upper Q 2 limit for the inelastic contribution, one should keep in mind 
that the basic representation with photon flux factorization is valid at Q 2 <C s only. Near 
this bound the original integrand is less than that obtained using our factorization, but its 
contribution to the total cross section is small. Therefore, the integration region can be 
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restricted from above by Q 2 <C s = yiV2S with sufficient accuracy. Since the virtual photon 
energies in both fluxes are roughly equal, we use 



01 



(18) 



as upper limit for the inelastic contribution without further reducing the accuracy. With 
this choice the factorization remains valid. 

In the y integration of ( |10"1) y max = 1 can be used as upper bound. The inaccuracy in the 
quantities that enter (structure functions) is inessential since contributions at yi ~ 1 can be 
safely neglected. 

III. TOTAL CROSS SECTION AND MASS LIMITS 

Elastic contribution. The proton form factors are written in the standard dipole approx- 
imation 



Hp 



G E (Q 2 



1 



[i + Q 2 /Q 



H P = 2.79285 , Ql = 0.71 GeV 2 . 



(19) 



Using dimensionless variables a = Am 2 /Ql and z = (a/4) xy 2 / (1 — y) the energy distribution 
is found[] as follows 



dN el _ y 2 {l-y) 
dy a 



a + z(l + fxl + Aa)) l(z,0) 



+ (a + z){ji 2 - l)I(z, a) - 



1 + z) 



(20) 



where 



I(z, a) 



L a + z * 1 fl - a 



With the numerical values for /i p and Qo the remaining y integration leads to the energy 
independent constant N e \ = 0.017672 as elastic contribution to N(E). 

Inelastic contribution. The structure functions F± t 2 are parametrized using the next to 
leading order quark distributions of [0 and the parton model relation F2 = 2xF\. To test 
the sensitivity of our results on this particular parametrization, we have repeated some of 



the calculations using the structure functions of fll6 |. The difference is small enough. 



The Q 2 dependence of the Fi j2 parametrizations is valid above some low input value Q 2 ow - 
Fortunately, the inelastic contributions to N(E) should be almost insensitive to the behavior 
of structure functions at small Q 2 . To test this statement, we consider two extrapolations 
for these functions below Qf ow : F 1:2 (x,Q 2 ) = or F h2 (x,Q 2 ) = F h2 (x,Ql w ). At E = 0.9 
TeV in GRV parametrization the results coincide within 1% accuracy. 



Compare fllil , Eq. (D.7) with a sign misprint corrected and upper Q 2 limit neglected. 
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Numerical estimates. At 0.9 TeV we obtain N(E) = 0.0410 in the GRV and N(E) = 
0.0338 in the MRRS approximation [|16| for this factor of Eq. (||). The quantity N(E) 
depends only weakly on the proton energy E. The ratio N(E)/N(0.9 TeV) varies from 0.966 
at E = 0.5 TeV to 1.006 at 1 TeV and 1.102 at E = 7 TeV using GRV parametrization. 

For a monopole with Jm = 1/2 and M/n = 1 TeV and a proton energy of 1 TeV we 
obtain the total cross section 

o- pp _ m x = 150fb. (21) 

It is useful to rewrite this photon production cross section for different proton collider ener- 
gies and different kinds of monopoles 



. MP , n) , WSP ^(^)\l^)\. (22 , 



Let us consider a luminosity integral of 2 fb _1 and a beam energy of 0.9 TeV (Tevatron). 
If we assume 10 events to be sufficient to detect the discussed effect, the following mass 
limits can be reached (for different spins Jm) 

( 0.998 TeV, J M = 
M < n ® \ 1.42 TeV, J M = 1/2 . (23) 
[ 2.56 TeV, J M = 1 

Taking 100 fb _1 and E = 7 TeV (LHC) we obtain the following mass limits 

f 7.40 TeV , J M = 
M <n®\ 10.5 TeV , J M = 1/2 . (24) 
[ 19.0 TeV , J M = 1 

The obtained limiting quantities correspond to a 77 — > 77 subprocess cross section 
(calculated as R (uj) & with an estimate of (u) taken from (|26| ) below) which is roughly 500 
pb for the Tevatron and 10 pb for LHC independent on monopole spin and n. These values 
are much higher than the cross section of the main background process 77 — > 77 via a W 



boson and a t quark loop which is about 20-30 fb [13]. 



IV. ENERGY AND MOMENTUM DISTRIBUTIONS 

Energy distribution for virtual photons. The photon fluxes in Eq. ([]) decrease with 
increasing photon energies. On the other hand, the 77 subprocess cross section rapidly 
increases with s = 4u)iu>2. Therefore, the main contribution to the cross section is given by 
region of intermediate cjj. As already mentioned, the dependence of the subprocess cross 
section on Q 2 can be neglected since the characteristic values of virtuality Q 2 Cs< M 2 . 
Therefore, the energy distribution for virtual photons is given by (compare Eqs. (|9.1C)) 



(]lfJ i,X ' 2 RE'ylf{ yi )ylf{y 2 ). (25) 



dy x dy 2 Vvr 
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Fig. Q shows the energy distribution y 2 f(y) for photons arising from one proton (the 
distribution is identical for each photon) using the two structure function parametrizations. 
The virtual photon energy distribution varies only weakly with the energy E, this weak 
energy dependence manifests itself in the weak E dependence of N(E) mentioned above. At 
0.9 TeV the average energy of the colliding photons and their energy spread are 

(u) = 0.314 E , (Aw) = 0.149 E . (26) 

77 differential cross sections. As noticed at the beginning of Sec. [TT], the transverse 
motion of incident (virtual) photons can be neglected with reasonable accuracy. Therefore, 
the transverse momenta of the produced photons are balanced: ptz ~ — Pta = Pt- The 
4-momenta of these photons in the c.m.s. of the protons can be written in two forms 

Pza = (£3,4, ±Pt, 0, PL3,Li) = Mcosh %,4, ±1, 0, sinh 773,4) . (27) 

Here transverse and longitudinal momenta and rapidities of photons are introduced. Using 
these notations we have 

With the standard transformation 

d _ s^Sid 



dyidy 2 dcos9 PTdpLzdp LA dpr 

the integrand of Eq. (|7]) (after integrating over Q 2 ) can be considered as the distribution 
over the momenta of produced photons. Then the differential cross section of 77 production 
via monopole loop can be presented in the form 



2 \ 2 / 1 

Pt \ _ L 1 



ujxuji J \ cosh [(773 - 774) /2], 



Integrating over one transverse momentum and azimuthal angle the cross section is given 

by 

d 3 a ( a\ 2 9 „. , 9 „, . 5RE 4 , .„ . 

ylf{yi)ylf{y2)— — (30) 



77 toto/ transverse momentum distribution. The total transverse momentum of the 
produced photon pair = Pt3 + PT4 is equal to the sum of transverse momenta of virtual 
photons, k<r = qri + <lT2- The latter are related to the photon virtualities by 

q^ = (l-^)(Q 4 2 -Q, 2 mi n), 2 = 1,2 (31) 

where Q 2 min is the minimal value of Q 2 for energy fraction as given in (O). 

Since characteristically Q 2 s, the photon pair transverse momentum is typically much 
smaller than the transverse momenta of the produced photons: <C Pt3 ~ Pta- Therefore, 
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the distribution in pu, pt factorizes from that in kx- The latter distribution can be obtained 
by integration over the virtual photon fluxes only (changing the order of integration). The 
detailed form of this dependence would be an additional test for the origin of the discussed 
photons. The corresponding calculations are simple but cumbersome, and one can postpone 
them to the time when first events of discussed type are observed. 

However, even before performing the calculation we can conclude that this distribution 
is peaked near kx = 0. This is evident for the elastic contribution where the scale of the 
distribution is limited from above by that of the form factor. For the inelastic contribution 
the virtual photon flux distribution is wider, however the mean value of k\ is much lower 
than s. 
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